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ABSTRACT 

We present spatially resolved (^50 pc) imaging of molecular gas species in the central kilopar- 
sec of the nearby starburst galaxy NGC 253, based on observations taken with the Atacama Large 
Millimeter/submillimeter Array (ALMA). A total of 50 molecular lines are detected over a 13 GHz 
bandwidth imaged in the 3 mm band. Unambiguous identifications are assigned for 27 lines. Based 
on the measured high C0/C 17 0 isotopic line ratio (>350), we show that 12 CO(l — 0) has moder¬ 
ate optical depths. A comparison of the HCN and HCO + with their 13 C-substituted isotopologues 
shows that the HCN(l-O) and HCO + (1-0) lines have optical depths at least comparable to CO(1-0 ). 
h 13 cn/h 13 CO + (and H 13 CN/HN 13 C) line ratios provide tighter constraints on dense gas properties 
in this starburst. SiO has elevated abundances across the nucleus. HNCO has the most distinctive 
morphology of all the bright lines, with its global luminosity dominated by the outer parts of the 
central region. The dramatic variation seen in the HNCO/SiO line ratio suggests that some of the 
chemical signatures of shocked gas are being erased in the presence of dominating central radiation 
fields (traced by C2H and CN). High density molecular gas tracers (including HCN, HCO + , and CN) 
are detected at the base of the molecular outflow. We also detect hydrogen /? recombination lines that, 
like their a counterparts, show compact, centrally peaked morphologies, distinct from the molecular 
gas tracers. A number of sulfur based species are mapped (CS, SO, NS, C2S, H2CS and CH3SH) and 
have morphologies similar to SiO. 

Subject headings: galaxies: formation — galaxies: evolution — individual galaxy (NGC 253) — 
astrochemistry; ISM — radio lines 


1. INTRODUCTION 

Characterizing the physical properties of the molecular 
gas phase in galaxies, the key phase for star formation 
(review by Kennicutt & Evans 2012), is among the main 
drivers in the studies of nearby galaxies. The 12 CO line 
is by far the brightest millimeter-wave molecular line in 
the interstellar medium (ISM) and has therefore been the 
main tracer of the molecular medium in both low- and 
high-redshift studies (e.g., reviews by Bolatto, Leroy & 
Wolfire 2013 and Carilli & Walter 2013). Even though 
CO emission is a reasonable tracer of the morphology 
and mass associated with the molecular ISM, observa¬ 
tions of other tracer molecules, that have different crit¬ 
ical densities and excitation temperatures, provide key 
information to constrain the physical processes in the 
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ISM. Studying a whole suite of molecular line tracers 
sheds light on the chemical state of the ISM, including 
the gas cooling and ionization balance in the molecular 
ISM. Ultimately, surveying provides the means by which 
to characterize molecular cloud conditions that are af¬ 
fected by the galactic environment, including feedback 
and shocks caused by star formation, and dynamical pro¬ 
cesses within a galaxy (e.g., Meier & Turner 2005). 

Millimeter line surveys of nearby galaxies using single¬ 
dish telescopes have shown that their spectra are rich in 
molecular lines (e.g., Usero et al. 2004, Martin et al. 2006, 
Costagliola et al. 2011, Snell et al. 2011, Aladro et al. 

2011, 2013). However, given the resolution of single-dish 
telescopes, most of these studies could only provide in¬ 
tegrated measurements, making it difficult to investigate 
spatial changes in chemical properties. Spatially resolved 
surveys using millimeter interferometers have been lim¬ 
ited to only a handful of transitions due to small band- 
widths and low sensitivity (e.g., Meier & Turner 2005, 

2012, Meier et al. 2014, Martin et al. 2014). These 
studies demonstrate the presence of strong chemical dif¬ 
ferentiation can exist within the molecular gas in nearby 
galaxies and their usefulness for constraining the evolu¬ 
tionary properties of the nucleus. 

The newly commissioned ALMA facility is revolution¬ 
izing chemical studies of galaxies. Its sensitivity, even 
in ‘Early Science’ mode is unprecedented, and molecules 
that are typically significantly fainter than the brightest 
ones (i.e., CO, HCN, HCO + , CN) can be detected in rea¬ 
sonable integration times. In addition ALMA’s large in¬ 
stantaneous bandwidth of Au =8 GHz implies that mul¬ 
tiple lines (in particular in the 3 mm band where Au/z/ is 
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highest) are covered in each observation. The study pre¬ 
sented here does not intend to discuss the details of every 
detected line but instead present the data and highlight the 
rich chemistry provided by ALMA’s new capabilities, in 
the nearby starburst galaxy NGC253. This paper com¬ 
plements two other studies of NGC253 using the same 
dataset, one describing the molecular outflow revealed in 
12 CO emission (Bolatto et al. 2013, hereafter B13) and 
the other constraining the molecular clump properties in 
the central starburst region using the brightest dense gas 
tracers available (Leroy et al. 2014, hereafter L14). 

The paper is structured as follows: Sec. [2] briefly sum¬ 
marizes the observations, data reduction and data prod¬ 
ucts. In Sec. [3] we present our identification of the rich 
suite of molecular line emission in the galaxy. We then 
proceed to sketch a simple picture of the conditions in the 
central region of NGC 253, by focusing on several groups 
of key molecular gas tracers (Se c J _[4|). We conclude by 
providing a short summary in Sec.JspThroughout the pa¬ 
per we adopt a distance to NGC 253 of 3.5 Mpc (Rekola 
et al. 2005), i.e., 1” corresponds to 17 pc at that distance. 
The systemic velocity of NGC 253 is ^250 km s -1 (e.g., 
Houghton et al. 1997). 

2. OBSERVATIONS AND DATA REDUCTION 

Observations: We have mapped NGC 253 with ALMA 
in cycle 0 (16 antenna array) in two different frequency 
setups and configurations. Details regarding the obser¬ 
vations and data reduction can be found in B13 and 
L14 and we here summarize the observational parameters 
that are most relevant for the current study. Both se¬ 
tups covered an instantaneous bandwidth of 8 GHz each: 
The first frequency setup, covering 85.6-89.6 GHz (lower 
side-band, LSB) and 97.4-101.4 GHz (upper side-band, 
USB), is a central 3-point mosaic along the major axis 
in the extended configuration (average beam size: ~2”, 
^35pc). This mosaic covers the central 1' (lkpc) of 
NGC253’s starburst. The second setup covers 99.8- 
103.7 GHz (LSB) and 111.8-115.7 GHz (USB) and is a 7- 
point mosaic of NGC 253’s center in the compact configu¬ 
ration. This yields an average beam size of ^4” (~70pc) 
and a field of view of roughly 1.5' (^1.5 kpc). 

Data reduction: All of the calibration and imaging of 
the data cubes was completed within CASA, including 
self-calibration to remove residual phase and flux calibra¬ 
tion errors. After imaging, much of the remaining analy¬ 
sis was performed with IDL scripts with some use of the 
Miriad and Matlab software packages. Most of the data 
presented here is based on the continuum-subtracted 
cube, where the continuum was defined in those spec¬ 
tral regions where no line emission was evident. Our 
final cube has 20kms _1 velocity resolution and a typical 
RMS of ^2.0 mJy beam -1 (extended, LSB: 1.9 mJy, ex¬ 
tended, USB: 2.1 mJy, compact, LSB: 2.0 mJy, compact, 
USB: 3.3mJy). 

Final data cubes: We display the information in the 
data cubes in two complementary fashions: through a 
position-velocity (pV) diagram along the major axis of 
NGC 253 in Fig. [l] a representation similar to long-slit 
spectroscopy in optical/NIR astronomy (Fig. [Th . The 
second (Fig. [2| is a single spectrum taken toward a cen¬ 
tral molecular peak (position 5). As outlined above, the 
data consist of four final data cubes, the LSB and USB 
cubes for the extended configuration, and those for the 


compact configuration. Figs. Hand [2] also show our line 
identification (see also Tab. m), as further discussed in 
Sec. OH 

Integrated line maps: We created line cubes for each 
line, and then blanked those with a mask derived from 
the CO(1-0) line (by far the brightest line in the band¬ 
pass). The integrated intensity map for a given line was 
then derived by simply adding all data in the respective 
data cube (without any further flux cutoff). Integrated 
line maps are corrected for primary beam attenuation 
and are shown in Fig. [3] 

Uncertainties: We calculated error maps by taking the 
number of channels per pixel into account. An inspec¬ 
tion of the data cubes revealed that these error maps 
provide too optimistic uncertainties, as they do not ac¬ 
count for artifacts (in particular the negative ‘bowl’ due 
to the missing short spacings) in the current data. We 
thus adopt a conservative 10% error for pixels that have 
been detected at high S/N (> 10), where the noise N 
is from taken from our error maps. For pixels that are 
detected at lower S/N (5<S/N<10) we adopt an even 
more conservative 30% uncertainty, to also account for 
possible issues in baseline determination. We ignore all 
pixels that have a S/N <5 in the analysis that follows. 
We present our line intensity measurements towards 10 
positions in NGC 253 (Fig. [4j coordinates in Tab. [TJ in 
Tab. [lj Toward the center of the galaxy (in particu¬ 
lar region 6, Fig. 0) line emission is observed against 
strong continuum (L14). Some resulting absorption will 
decrease the integrated line signal at these locations. Be¬ 
cause of this the fluxes toward this location should be 
considered highly uncertain. 

Spatial filtering: As ALMA is an interferometer it acts 
as a spatial filter, sampling only a range of spatial scales. 
So the observations will potentially resolve out some flux. 
In the compact (extended) configuration spatial scales of 
12 - 47 kA (23 - 100 kA) [90th percentile] were sampled, 
corresponding to 4.4" - 18" (2" - 9"). Hence the observa¬ 
tions should adequately sample fluxes uniform over ^10" 
(in one channel). The one exception to this is 12 CO(1-0), 
which has been zero-spacing corrected and therefore de¬ 
tects all flux (see discussion in B13). It is not possible to 
estimate the percentage of detected flux for every line in 
the survey because suitable single-dish observations are 
often not available. We do, however, determine this frac¬ 
tion for a number of transitions where possible. For the 
following transitions we calculate detected flux percent¬ 
ages over single-dish beam of ^22 " - ~28" of: HCN(1- 
0) — 100 % (Paglione et al 1995), HCO+(1-0) — 55 % 
(Martin et al. 2009), C 17 O(l-0) — 60 % (Henkel et al. 
2014), SiO(2-l) — 110 % (Martinet al. 2009), CN(l-0;§ 
- \) — 60 % (Henkel et al. 2014), C 2 H(1-0;| - \) — 60 
% (Nakajima et al. 2011) and HNCO(4 0 4-3o3) — 50 % 
(Nguyen-Q-Rieu et al. 1991) [though this value is highly 
uncertain because HNCO does not peak where Nguyen- 
Q-Rieu et al. pointed]. Therefore, it appears that the 
data consistently detect at least 50 % of their respective 
single dish fluxes. Furthermore our discussion focuses on 
the compact clumps of emission, where much higher frac¬ 
tions of the flux are detected (->90 %). (In fact, spatial 
filtering of the interferometer actually affords advantages 
because it allows the separation of these compact struc¬ 
tures from any extended diffuse medium that can domi- 
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TABLE 1 

Observed Emission Lines 


Transition 

Rest Frequency 
(GHz) 

Setup 

Comments 

HC 15 N(J = 1-0) 

86.054 

ext, LSB 

tent. ID 

SO(J jv =2 2 -1i) 

86.09395 

ext, LSB 


H 13 CN(l-0) 

86.34016 

ext, LSB 


HCO(Jiy a iy c =li,o-Oo,o) 

86.67076 

ext, LSB 


H 13 CO+( 1 - 0 ) 

86.75429 

ext, LSB 


SiO(2-l; u=0) 

86.84696 

ext, LSB 


HN 13 C(l-0) 

87.09085 

ext, LSB 


C 2 H(JV=l-0;J=3/2-l/2) 

87.31690 

ext, LSB 

main C 2 H FS component 

C 2 H(l- 0 ;l/ 2 -l/ 2 ) 

87.40199 

ext, LSB 

secondary C 2 H FS component 

u 

87.567 

ext, LSB 

poss. ID: SiN( 2 - 1; 5/2 - 3/2) 

CH 3 NH 2 ( 3i -3o) 

87.785 

ext, LSB 

tent. ID 

NH 2 CHO( 4i, 3 -3i, 2 ) 

87.853 

ext, LSB 

tent. ID 

HNCO(4 0 ,4 - 3o, 3 ) 

87.92524 

ext, LSB 


H(52 )/3 

88.40569 

ext, LSB 


HCN(l-O) 

88.63160 

ext, LSB 


HCO+( 1 - 0 ) 

89.18852 

ext, LSB 


34 SO(3 2 - 2 j) 

97.720 

ext, USB 

tent. ID 

CS(2-1) 

97.98095 

ext, USB 


CH 3 CH 2 CN 

98.525 

ext, USB 

tent. ID; blend: (116 — IO 7 ) & (II 7 — IO 7 ) 

u 

98.667 

ext, USB 

no obvious ID 

H(40 )a 

99.02295 

ext, USB 


H(50 )/3 

99.22521 

ext, USB 


SO(3 2 — 2i) 

99.29987 

ext, USB 


U 

99.669 

ext, USB 

no obvious ID 

C 2 S( 7s - 67 ) 

99.86651 

ext, USB 


HC 3 N( 11 - 10 ) 

100.07640 

ext, USB, comp, LSB 


u 

100.542 

ext, USB, comp, LSB 

poss. ID: CH 2 CN( 5 2 3 - 4 2 2 ; 9/2 - 7/2) 

u 

100.612 

ext, USB, comp, LSB 

poss. IDs: CH 2 CN( 5 0 5 - 4 0 4 ; 9/2 - 7/2) 




CH 3 CH 2 CN(lh j 0 - 10i, 9 ) 

u 

100.632 

ext, USB, comp, LSB 

poss. IDs: NH2CN(5 i 4 — 4i 3 ); 




CH 2 CN{5 2 4 - 4 2 ’ 3 ; 11/2 - 9/2) 

u 

100.992 

ext, USB, comp, LSB 

poss. ID: CH 3 CH 2 OH( 82,7 - 8is) 

CH 3 SH(4o - 3o)A+E 

101.13916 

ext, USB; comp, LSB 


CH 3 SH( 4 2 - 3 2 )A 

101.180 

ext, USB; comp, LSB 

tent. ID 

H 2 CS(3i, 3 - 2i, 2 ) 

101.47788 

comp, LSB 

ortho-Lucs' 

CH 2 CO( 5i, 4 -4 4i3 ) 

101.892 

comp, LSB 

tent. ID 

U 

101.988 

comp, LSB 

no obvious ID 

NH 2 CHO(5i,5 — 4i >4 ) 

102.070 

comp, LSB 

tent. ID 

CH 3 C 2 H(6 fc - 5 fc ) 

102.54798 

comp, LSB 


h 2 cs 

103.04054 

comp, LSB 

para: (32,2 — 22 , 1 ), (3o,3 — 20 , 2 ), (32,1 — 22 , 0 ) 

C 2 S(88 — 77 ) 

103.64075 

comp, LSB 


CH 3 CHO(6 1 , 5 -5 1 , 5 )A+ 

112.247 

comp, USB 

tent. ID 

C 17 O(l-0) 

112.35928 

comp, USB 


HCOOH 

112.459 

comp, USB 

tent. ID, blend: ( 53,3 — 43 , 2 ) & (53,2 — 43,i) 

CN(l-0;l/2-l/2) 

113.19128 

comp, USB 

secondary CN FS group 

CN(l-0;3/2-l/2) 

113.49097 

comp, USB 

main CN FS group 

U 

114.218 

comp, USB 

poss. IDs: CO(1-0; v=l); 




C 4 H(12-11; 25/2-23/2& 23/2-21/2) 

U 

114.605 

comp, USB 

no obvious ID 

CH 3 CHO( 6 o, e ~ 5 0 , 3 )A++E 

114.952 

comp, USB 


NS(5/2-3/2;7/2-5/2) 

115.15394 

comp, USB 

blended with CO 

CO (1-0) 

115.27120 

comp, USB 

Bolatto et al. 2013 

NS(5/2-3/2;7/2-5/2) 

115.55625 

comp, USB 



Note. — Note: species in italics represent tentative identifications. In these cases, the frequency listed is the observed one corrected to rest 
frequency based on the known velocity field and does not reflect the rest frequency of the tentatively identified line. Lines labeled U are the ones 
that we could not securely identify, and in some cases we give possible ID’s in the last column. 


nate single-dish observations.) Finally, we avoid compar¬ 
ing line intensities derived from different array configura¬ 
tions wherever possible in this study, to further mitigate 
against differences in resolved flux. Hence we conclude 
that uncertainties in line ratios due to different degrees 
of missing flux are <10 %. 

3. RESULTS 

3.1. The Nucleus of NGC253 


The nucleus of NGC 253 is characterized by the inner 
portion of NGC 253’s large-scale bar: The highly inclined 
(i ~ 78 deg) nuclear disk extends from the very center out 
to a radius, r > 370 pc. In some transitions emission is 
seen out to the edge of the mosaic. The outer part of 
the nuclear disk (the ‘outer nuclear disk’) is suggested to 
represent the location where gas flowing radially inward 
along the large-scale bar collects between the outer and 
inner Lindblad resonance (Garcia-Burillo et al. 2000). 
Embedded within the outer disk is a compact (r ~ 170 
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Fig. 1.— Position-velocity cut along the major axis of the central bar of NGC253 (equivalent to an optical/near-infrared long-slit 
spectrum). Angular offsets are given with respect to the galaxy center. The left two panels show the lower and upper sidebands of the 
extended ALMA configuration, the right two panels the corresponding sidebands of the compact configuration. Identified molecular spectral 
lines are labeled (secure: bold, tentative: italics, see Tab. ^for a full description of transitions). The CO line is completely saturated in this 
presentation and the HCN, HCO+ and CN lines are partially saturated. Table ^ lists the coordinate about which the offsets are referenced. 


TABLE 2 

NGC 253 Positions 


Region 

RA(2000.0) 

DEC(2000.0) 

1 

00 47 33.041 

-25 17 26.61 

2 

00 47 32.290 

-25 17 19.10 

3 

00 47 31.936 

-25 17 29.10 

4 

00 47 32.792 

-25 17 21.10 

5 

00 47 32.969 

-25 17 19.50 

6 

00 47 33.159 

-25 17 17.41 

7 

00 47 33.323 

-25 17 15.50 

8 

00 47 33.647 

-25 17 13.10 

9 

00 47 33.942 

-25 17 11.10 

10 

00 47 34.148 

-25 17 12.30 

Map 

00 47 33.100 

-25 17 17.50 


Note. — Coordinates of the 10 regions in NGC 253 (indicated in 
Fig. pj for which intensities have been measured (see Tab. ^). The 
coordinates of (0, 0) reference position of the maps is also given. 


pc) region exhibiting a large quantity of high density gas 
and intense star formation, which is discussed in detail 


in Section 4.2 A molecular gas outflow/wind is being 


driven from this inner nuclear disk, as discussed in B13. 

We have selected 10 individual locations that span the 
nuclear disk and the base of the molecular outflow in 
NGC 253 (Tab. [l]) . Regions 3-9 trace the inner starburst 
disk from west to east and regions 5 and 6 are coincident 
with the detected millimeter continuum emission in the 
galaxy (L14). Regions 1 and 10 are located at the west¬ 
ern and eastern base of the southern molecular outflow 
(B13). Region 2 indicates the shocked region towards the 
north. The GMC physical properties in these regions are 
discussed in detail in L14. We further caution that some 
of our line measurements towards region 6 (the center) 
may be affected by absorption and continuum subtrac¬ 
tion uncertainties (as discussed in Sec. [2]). 

In Tab. [T| we provide a summary of all detected lines 
(column 1)7 their rest frequencies (column 2) and obser¬ 
vational setup (column 3). For each of our 10 regions, 
we measure the peak integrated intensity (in Kkms -1 ) 
at that position and report the values in Tab. [l] We 
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Fig. 2. — Extracted spectrum near the center of NGC253 (position 5 in Fig. [ 4 J, on a logarithmic scale. To deal with subtle issues 
regarding the bandpass calibration, the data have been independently continuum subtracted for each spectral window. After that, the 
average continuum flux density in a given spectral window was added to the spectrum. The flux density offset between the observations 
obtained in the compact and extended configurations is real; the observations in the compact configuration recover more of the extended 
flux. As in Fig.^ identified molecular spectral lines are labeled (secure: bold, tentative: italics, see Tab. ^for full description of transitions). 


measure all values from beam-matched maps, i.e., all the 
compact and extended configuration data have been con¬ 
volved to a common resolution, respectively (this com¬ 
mon resolution is given in column 4 of Tab. [l] whereas 
column 3 gives the original resolution). We have dis¬ 
cussed our uncertainty estimates in Sec. [2] Tab. [l] also 
includes the intensity measurement in each transition for 
the entire galaxy (columns 5 and 6). 

3.2. Molecular Emission Across NGC 253’s Nucleus 

Fig. [l] shows the position-velocity (pV) diagram ori¬ 
ented along the molecular bar (i.e., along the major axis 
of the galaxy) of all four data cubes. This figure shows 
the richness of the ALMA data cubes, and the challenges 
in finding the correct line identification. Already from 
this presentation it is clear that the emission of the dif¬ 
ferent line species emerge from different regions in the 
galaxy. 

We detect at least 50 emission lines in total of which 
27 have unambiguous identifications. This corresponds 
to a detection ratio of one line every ^0.25 GHz. Species 
printed in italics are tentative IDs (13 lines). These lines 
have clear candidate transitions from a species no t ye t 
well established in extragalactic systems (see Sec. _|3.3| ). 
The candidate IDs are given in column 4 of Tab. [l] In 
some cases we were not able to find a plausible line identi¬ 
fication or the broadness of the spectral line prohibited a 
unique identification. Those lines are labeled U (uniden¬ 
tified; 10 cases). In those cases we give our best estimate 
for each feature’s rest frequency in column 2. 

We show line maps for 27 of the bright, unblended lines 
in Fig. [3j In Fig. [4] we select out six species that summa¬ 
rize the basic morphological patterns seen in the sample. 
Each represent key tracers of the different phases of the 
molecular gas, as discussed below. These are: 12 CO, a 
tracer for the overall molecular gas distribution, 12 C 17 0, 


an optically thin tracer of molecular gas emission, C 2 H, a 
photon-dominated region (PDR) tracer, H40cq a hydro¬ 
gen recombination line and H II region tracer, HCN, a 
high-density gas tracer and HNCO, a weak shock tracer. 

From Figs. [3] and [4] it is apparent that the morpholo¬ 
gies change significantly between species (as is already 
evident from Fig. [I]). Molecular gas column densities re¬ 
main quite large across the entire nuclear disk, maintain¬ 
ing N(H 2 ) > 10 23 cm -2 over much of the disk (Sec. 3.4). 
The wide range of dense gas tracers observed demon¬ 
strates that dense gas is present across much of this re¬ 
gion, but the dense gas fraction increases toward the in¬ 
ner disk (L14). The millimeter hydrogen recombination 
lines, including for the first time H/3 lines, show that 
dense H II regions associated with the young starburst 
(or possibly AGN; Mohan et al. 2002) are confined to 
the innermost part of this inner disk. Species like C 2 H, 
CN and CHsC 2 H also dominate from this inner nuclear 
disk. HNCO, on the other hand, is dominated by the 
outer nuclear disk. 

The brightest and most widespread transitions (HCN, 
HCO + and CN) show weak emission extended vertically 
between 1 and 2, following the base of the molecular 
outflow, suggesting there is dense molecular gas here. 
However, we do not discuss the dense gas tracers in the 
outflow in detail here due to the difficulties dealing with 
the morphology of the maps above and below the disk in 
the presence of the negative ’bowl’. 


3.3. Tentative and Unidentified Lines 

There are a number of detected lines for which no clear 
identification was possible. We have listed tentative iden¬ 
tifications in Tab. [I]- in some cases no plausible species 
were found, and these lines constitute unidentified ( U) 
lines. The majority of these tentative identifications (but 
not tentative detections) match prominent transitions 
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Fig. 3.— Integrated intensity maps of selected molecules in our observations (Tab. IT]). The respective intensity scale (in units of Kkms -1 ) 
is shown in each panel. The name of the line and its rest frequency are given in the top right corner of each panel. The beam size is shown 
in the bottom right corner. Tentatively identified lines are labeled in italics. 

from larger molecules, including the aldehydes CH3CHO 
and NH2CHO, cyanides CH 2 CN and CH 3 CH 2 CN and 
a number of organics of similar complexity, CH 2 CO, 

HCOOH and CH 3 CH 2 OH. These molecules are abun¬ 
dant in the Galactic center (e.g., Cummins et al. 1986) 
and are expected to be detectable at our sensitivity in 
NGC253. If follow-up, multi-line studies confirm these 
IDs then they would represent their first extragalactic 
detections. 

3.4. Excitation and Abundances 

When appropriate, an estimate of abundances is made 
for the detected species. Molecular column densities 
Nmol are determined assuming optically thin LTE emis¬ 
sion: 


where Sul the line strength, and g u 's and E u are the 
upper state degeneracy and energy, respectively, Q is 
the partition function, /i 0 is the dipole moment in debye 
and T ex is the excitation temperature associated with 
the transition (e.g., Turner 1991). For the symmetric 
(CH 3 SH) and asymmetric (HNCO, H 2 CS) tops Q is pro¬ 
portional to Te^ 2 , whereas Q for the linear rotors (the 
rest) is proportional to T ex . 

To calculate abundances it is necessary to adopt an 
T ex . Since this survey covers only one band, gas exci¬ 
tation is not the focus of this project. The only lines 
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that are detected and directly constrain gas excitation is 
a pair of SO transitions, SO(32 — 2i) and SO(2 2 — li). 
From this pair, an excitation temperature is estimated, 


Tex = 


10.09 K 


( 2 ) 


where R is the intensity ratio, SO(3 2 — 2 i)/SO( 2 2 — li). 
Towards the inner nuclear disk this ratio has a value of 
R~2.5. There is some evidence for the ratio to increase 
towards the outer nuclear disk (R<4), but the faintness 
of SO(2 2 — li) limits what can be said. For these ratios, 
we estimate the inner disk to have an excitation temper¬ 
ature of 74 K in SO, dropping to T ex <17 K. Indirect 
methods of estimating T ex can also be a chieved by ra¬ 
diative transfer modeling (see section 4.2). There T ex is 
found to range between ^6 - 45 K. Therefore, T ex likely 


falls between at least 10 - 75 K, across much of the nu¬ 
clear region. Ott et al. (2005) use NH3 to estimate gas 
kinetic temperatures, T^ n , across the inner nuclear disk 
of approximately twice the 74 K measured here, however 
Knudsen et al. (2007) (and this paper) favor T k in nearer 
this value for HCN and HCO + (see section 4.2). The 
drop in T ex toward the outer nuclear disk region reflects 
SO becoming more strongly sub-thermal in lower density 
gas. 

The relevant excitation temperature for a given species 
/ transition depends on a number of factors, including 
gas density, opacity, molecular structure and location 
within the nucleus. For transitions with effective critical 
densities significantly above the gas density, T ex will be 
lower than the T kin (subthermal excitation). Changes in 
tih 2 and T^ n with position in the nucleus imply changes 
in excitation. The effective critical density of a species 
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depends on line opacity, r, (photons become an effec¬ 
tive ’collision’ partner) and hence abundances, with high 
opacity lowering critical densities, n cr , roughly in pro¬ 
portion to r. Finally, even in the case that a transition 
is thermalized, the LTE calculated columns are strict up¬ 
per limits because it assumes that all transitions of the 
molecule are thermalized up to arbitrarily high energies. 

With these caveats stated, abundances are the impor¬ 
tant astrochemical variable and so we report ’reference’ 
ranges. In Table [lj abundances are reported for T ea , 
ranging from 10 K (first number in each entry) to 75 K 
(second number). Also included in the the table are n cr 
(not including opacity effects) for the transitions. It is 
expected that for the inner nuclear disk positions tran¬ 
sitions with moderate n cr (e.g., C 17 0, HNCO, H2CS, 
C2H, HCO + , and SO) will have T ex at the high end of 
the range. High n cr transitions and outer nuclear disk 


locations likely will have T ex (and hence abundances) at 
the low end of the range. Future, follow up, multi-line 
studies of different species are necessary to narrow these 
quoted ’reference’ value ranges. 

To convert N mo i to fractional abundances we also re¬ 
quire H2 column density, N(H2). N(H2) is most easily 
obtained from the CO(1-0) brightness and an empirical 
conversion factor, Xqo- We adopt a CO-to-H 2 conver¬ 
sion factor of Xco = 0.5 x 10 20 cm -2 (K km s -1 ) -1 . 
This is a factor of four lower than typical for the Galac¬ 
tic disk (Strong et al. 1988; Hunter et al. 1997; Bolatto 
et al. 2013a), but consistent with what has been previ¬ 
ously estimated for NGC 253 (e.g. Paglione et al. 2001, 
B13, LI4). Uncertainties in Xc o ar e significant, likely 
at the ±2x level, and in section |4T| some discussion of 
its validity in the context of optical thin C 17 O(l-0) is 
discussed. 







ALMA Multi-line Imaging of the Nearby Starburst Galaxy NGC253 


9 



30 20 1 0 0 -10 -20 -3030 20 1 0 0 -10 -20 -3030 20 1 0 0 -10 -20 -30 

R.A. Offset (") R.A. Offset (") R.A. Offset (") 


Fig. 4.— Collection of key molecular gas tracers in NGC253. From top left to bottom right: CO(l-O) (main molecular gas tracer), 
C 17 O(l-0) (optically thin tracer of the molecular gas), C 2 H (PDR tracer), H40ai (radio recombination line), HCN (high density gas tracer), 
HNCO (shock tracer). The crosses and labels (bottom left panel) mark the different regions defined in the galaxy. Logarithmic contours 
start at lKkms -1 and increase/decrease (white/grey) by a factor of 3 from one contour to the next. The beam size is indicated in the 
bottom left corner of each panel. 


4. DISCUSSION 

In the following we present a broad discussion of the 
chemistry in the nuclear region of NGC 253. The primary 
goal is to obtain spatially resolved views of the different 
molecular environments in and around the starburst cen¬ 
ter. 

We begin our discussion by studying the main tracer of 
the molecular gas, 12 CO, its isotopologues, and implied 
molecular gas opacity in Sec. |4 . 11 In Sec. 42] we con¬ 
tinue with a similar discussion of the isotopologues of the 
main tracers of the dense molecular gas phase (H 13 CN, 
H 13 CO + , HN 13 C). Thi s is followed by a discussion of 
tracers of PDRs, (Sec. 4.3), shock tracers (HNCO and 
SiO; Sec 
bination 


4.4| ), Sulf ur sp ecies (Sec. 4.5), hydrogen recom- 
mes (Sec. 


4.6| ), and other tentatively identified 
species (Sec. |4.7| ) detected in our study. A schematic 
picture summarizes the overall chemistry of the central 
region of NGC 253 in the conclusions/summary section 
(Sec.§. 


4.1. CO Isotopologues and Gas Opacity 

Our data cover a number of isotopologues (i.e. 
molecules where one atom is replaced with an isotope). 
Given the fact that atomic isotopes are much less abun¬ 
dant than the main atom, many of these isotopologue 
lines have low optical depth compared to the main 


species. Isotopologue ratios thus provide insights into the 
optical depth of the main lines and can also be used to 
constrain isotopic abundance ratios. The isotopic abun¬ 
dance ratio ultimately can provide insights into stellar 
nucleosynthesis and possible variations of the stellar ini¬ 
tial mass function (IMF). In the discussion that follows 
we do not interpret the sometimes anomalous line ratios 
found for region 6 (i.e. the region that is coincident with 
the central continuum emission) as it is likely that in 
some cases our measurements are affected by absorption 
(Sec. [2). 

A comparison of the brightness of the C 17 O(l-0) tran¬ 
sition to the 12 CO(1-0) transition [hereafter CO(1-0)] 
can provide important constraints on the total H2 col¬ 
umn. Being an optically thin version of CO, even in 
this extreme star formation environment, it permits CO 
abundances to be determined by ‘counting molecules’ if 
the [CO/C 17 0] abundance ratio is known. 

The [CO/C 17 0] abundance ratio is constrained based 
on existing measurements and nucleosynthetic expecta¬ 
tions. In the Galaxy, [CO/C 17 0] is <1900 at the solar 
radius and drops to ^900 in the Galactic center (e.g., 
Wilson & Rood 1994, Ladd 2004, Wouterloot et al. 2008). 
This decrease in the [CO/C 17 0] abundance is consistent 
with 17 0 being a secondary chemical evolution product 
formed from reactions between primary 16 O and a pro- 
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ton in intermediate mass stars. Therefore, in a strongly 
processed location like the center of a starburst, a low 
value of the C0/C 17 0 abundance ratio is expected, so 
we adopt [C0/C 17 0] = 1000, which is approximately the 
Galactic center value referenced above, as our nominal 
value. 

As shown in Fig. [5| this isotopologue line ratio is fairly 
constant with position in NGC 253 and large (>350). For 
a ratio 16 0/ 17 0 of 1000 in NGC 253, the measured ratio 
nominally implies a 12 CO opacity of ~2.5, a moderate 
value for the main CO line. Raising the intrinsic 16 0/ 17 0 
abundance to 2000 gives a 12 CO opacity of order ^5.5. 

Various non-LTE effects can occur in a starburst en¬ 
vironment that could alter the observed CO/C 17 0 line 
ratio. Examples that artificially inflate the CO/C ir O 
line ratio include: 1) CO has a higher excitation tem¬ 
perature than C 17 0 because CO has high opacity and 
its emission is dominated by warmer, externally heated 
edges of the GMCs while C 17 0 is dominated by the 
cooler interiors. 2) Isotope-selective photo-dissociation 
makes the cloud sizes smaller in C 17 0 than CO. Or, 
3) due to the absence of radiative trapping in C 17 0, 
C 17 0 is subthermal relative to CO and therefore fainter 
(e.g., Meier & Turner 2001). The first two would be 
expected to be pronounced close to the nuclear star for¬ 
mation. However, interestingly, no evidence is seen for 
an increase in the CO/C 17 0 rati o to wards the inner disk, 
where PDRs are enhanced (Sec. 4.3) and warm gas dom¬ 
inates. But if there is a decrease m density away from 
the inner disk then the third option can operate in the 
outer disk, partially compensating for the first two. For 
CO/ 13 CO abundances [CO/ 13 CO] > 60 (Martin et al. 

2010) , CO opacities estimated from the CO/C 13 0 line 
ratios are ^5 (e.g., Paglione et al. 2004, Sakamoto et al. 

2011) . This scenario is consistent with CO/C 17 0 data 
(and CO/C 18 0) if non-LTE effects are at most a factor of 
~2. We conclude that, despite the large column density 
towards NGC 253’s nucleus, CO(l-O) is only modestly 
opaque. It is likely that this low optical depth is due to 
increased line widths due to turbulence, the molecular 
outflow, non-circular motions and/or elevated gas tem¬ 
peratures that reduce the CO opacity per unit velocity. 

The moderate CO optical depth is also reflected in the 
implied sub-Galactic CO-to-H 2 conversion factor in star- 
bursts like NGC 253 (Harrison et al. 1999, Downes & 
Solomon 1998, Meier & Turner 2001; Bolatto, Leroy & 
Wolfire 2013). A comparison of the LTE C 17 0 column to 
the H 2 column estimated by an Xco of 0-5 x 10 20 cm -2 
(K km s -1 ) -1 , gives a value for [C 17 0/H2] constant (to 
within 20%) of 1.9 x 10 -7 across positions 3-10 (using 
Eq. 1 and T ex =75K, Sec. [IZ]). If [CO/C 17 0] = 1000 
and [CO/H 2 ] = 10 -4 , then we would expect a value of 
[C 17 0/H2] — 10 -7 . The level of agreement between the 
two methods of constraining N(H 2 ) is reasonable, being 
less than the uncertainties in each method. The small 
differences are an indicator of the uncertainty in 7V(iJ 2 ). 
Agreement could be made exact by raising [C 17 0/H2] 
and/or by lowering the adopted T ex . 


4.2. Dense Gas Tracers and their Isotopologues 

The commonly observed dense gas line ratios, 
HCO+/HCN, HNC/HCN and HCO+/HNC ratios are 


predicted to be sensitive to gas density, cosmic ray ion¬ 
ization and possibly the X-ray versus UV ionization rate 
or the present of mechanical heating (e.g., Kohno et 
al. 2001, Meijerink & Spaans 2005, Gracia-Carpio et al. 
2006, Knudsen et al. 2007, Papadopolous 2007, Krips et 
al. 2008, Baan et al. 2008, Loenen et al. 2008, Meijerink 
et al. 2011, Kazandjian et al. 2012). HCO + (1-0) ef¬ 
fective critical densities are nearly a factor of five lower 
than HCN, so to first order, elevated HCO+/HCN and 
HCO+/HNC ratios are expected in moderate density gas 
(' n H 2 10 4,5 cm 3 ). As densities increase, both transi¬ 
tions thermalize and the ratio tends to unity. Moreover, 
HNC/HCN ratios are expected to drop in energetic en¬ 
vironments because HNC is preferential destroyed in hot 
gas (e.g., Schilke et al. 1992, Meier & Turner 2005). If 
these transitions have high optical depths, which is prob¬ 
able given the large columns of dense gas, all ratios are 
driven toward unity. As a result of the large inferred 
optical depths, changes in these main isotopic line ra¬ 
tios due to interesting physical/chemical changes will be 
partially hidden. 

Here we focus on the dense gas properties in the cen¬ 
tral starburst region of NGC 253. The HCO+/HCN 
intensity ratio (Fig. [5| is observed to be slightly less 
than unity and exhibits very little change over the in¬ 
ner disk (except for position 2 and 3). This sug¬ 
gests that both HCN(l-O) and HCO + (1-0) are optically 
thick and slightly subthermal across much of the nucleus 
(Knudsen et al. 2007). Here we add new HCN/H 13 CN, 
HCO+/H 13 CO+, H 13 CO+/H 13 CN, and HN 13 C/H 13 CN 
intensity ratios to further constrain dense gas properties. 
The isotopologues of HCN, HCO + , and HNC have lower 
optical depths and allow estimates of the line opacity and 
changes in the abundance ratio of these species. 

The observed HCN/H 13 CN and HCO+/H 13 CO+ both 
have values of 10 - 15 over much of the inner nuclear disk 
(Fig# The fact that these ratios are much lower than 
the expected [ 12 C/ 13 C] ratio of >60 (Martin et al. 2010) 
clearly demonstrates that both HCN(l-O) and HCO + (l- 
0) have opacities greater than one. This confirms the 
conclusions of earlier work (Knudsen et al. 2007). For 
A = 12 C/ 13 C abundance ratios of 80 above and LTE as¬ 
sumptions, both HCN and HCO + have r ~ 5 — 8. Inter¬ 
estingly this is an optical depth similar or slightly larger 
than inferred for CO(l-O) from C 17 0. This is expected 
because the HCN/H 13 CN and HCO + /H 13 CO + ratios are 
quite similar to the observed CO/ 13 CO line ratios (e.g., 
Paglione et al. 2004; Sakamoto et al. 2011). Such optical 
depths for HCN and HCO + are large. For a T ea ,=75 K, 
inferred optically thin HCN abundances are ^ 2.5 x 10 -8 
(Table 4). However, if the isotopologues of HCN, HCO + 
and HNC have a lower excitation temperature than the 
main species (likely true, see below), then LTE inferred 
opacities are lower limits. 

The H 13 CO+/H 13 CN ratio in particular shows larger 
variation (Fig. [5| than seen in the main isotopic ratios. 
Moreover the average values of the ratios are different. 
Toward the two main inner disk GMCs (positions 5 and 
7) HCO+/HCN and H 13 CO+/H 13 CN have very similar 
ratios, but toward the outer nuclear disk the isotopic 
dense gas ratios decrease relative to their main versions. 
Likewise the average HN 13 C/H 13 CN ratio of ^0.35 is 
significantly lower than the single-dish value for the main 
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Fig. 5.— CO and dense gas line ratios (Sec s. |4.1|a nd j4.2|>. The respective line ratio is given in each panel. The x-axis number corresponds 
to the region number as defined in Sec. |3.2| and Fig. El The horizontal line in each panel shows the value for the entire galaxy (i.e., the 
global ratio). Uncertainties are determined from TablefD 


species ratio (~0.74; Baan et al. 2008). The fact that the 
isotopologue substituted species have ratios farther from 
unity is consistent with the expectation that high optical 
depth in the main lines artificially drives the ratio closer 
to unity. These isotopologue ratios tend to move NGC 
253 further into the PDR-like ratio parameter space of 
Baan et al. (2008), consistent with dense gas abundance 
ratios influenced by warm, PDR across the inner nuclear 
disk. 

Given the lower opacities of the isotopologues of HCN, 
HCO + , and HNC, it is expected that subthermal exci¬ 
tation will be the most relevant non-LTE consideration, 
therefore it is worth modeling the physical conditions 
implied by these lines. Knudsen et al. (2007) have mod¬ 
eled the main isotopologues, including selected AJ tran¬ 
sitions. Here we carry out an independent large velocity 
gradient (LVG) radiative transfer modeling focusing on 
the isotopic line ratios. Line ratios are modeled as they 
are to first order independent of the unknown areal fill¬ 
ing factor, and therefore constrain the parameter space 
more precisely. The model used is patterned after the 
models of Meier et al. (2008), with J max extended to 
20 and collisional coefficients adapted from the Leiden 
LAMDA database (van der Tak et al. 2007; see Table 
4 for references to the individual molecular rate coeffi¬ 


cients). To model the radiative transfer, assumptions 
must be made about the abundance per velocity gradi¬ 
ent. Models with a wide range of abundance per velocity 
gradient (10 -11,33 — 10 -7,67 km -1 s pc) were calculated. 
In general it is not easy to obtain agreement for the line 
ratios of HCN/HCO+, HCN/H 13 CN, HCO+/H 13 CO+, 
H 13 CO+/H 13 CN, HN 13 C/H 13 CO+ and HN 13 C/H 13 CN. 
HCO +, s effective critical density is significantly lower 
than HCN and so HCO + abundances have to be low com¬ 
pared to HCN to obtain ratios of HCO+/HCN lower than 
unity. But pushing HCO + abundances too low causes 
disagreement with the isotopic ratios. If we focus only 
on the five isotopic line ratios, then reasonable solutions 
that match all are found over a fairly narrow range of pa¬ 
rameter space. We adopt Xhcn /dv/dr = 10 -8,1 km -1 s 
pc and Xhco+/ dv/dr = Xhnc /dv/dr = 10 -8 - 85 km -1 
s pc, consistent with Xhcn — 2.3 x 10 —8 , Xhco+ — 
Xhnc — 4.2 x 10 -9 and a velocity gradient of ^3 km 
s -1 pc -1 . For the isotopologues, a 12 C/ 13 C ratio of 80 is 
adopted. These values are in reasonable agreement with 
the values listed in Table 4 and the observed gas kinemat¬ 
ics. There is room to adjust the abundances per veloc¬ 
ity gradient of the five species individually, but keeping 
HCO + and HNC abundances equal does an acceptable 
job. The low abundance of HCO + demanded relative to 
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Fig. 6.— LVG modeling of the inner nuclear disk. The acceptable model predicted parameter space for (clockwise) HCN/HCO+ 
HCN/H 13 CN, HCO+/H 13 CO+, H 13 CO+/H 13 CN, HN 13 C/H 13 CO+, and HN 13 C/H 13 CNl-0 transitions. Contours are the Tier observed 
line ratio for the average of position 5 and 7. Tick marks indicate the direction of the favored pa rame ter space if it is not obvious. Adopted 
abundance per velocity gradient values from the six different species are discussed in sectio n 14.2| The cross marks an example set of 
parameters (T/^ n = 80 K, tih 2 — 10 4-75 cm -3 ) that is within la of with all five isotopic ratios. The HCN/HCO+ line ratio is not matched, 
but is within 30 % of the predicted ratio. 


HCN is a direct consequence of its lower effective critical 
density. 

Figure [ 6 ] displays the acceptable parameter space for 
the ratios observed toward the inner nuclear disk (av¬ 
erage of the ratios toward positions 5 and 7) displayed 
in Figure [5] When admitting uncertainties, good agree¬ 
ment to all five ratios is obtained for n # 2 = io 4 - 5-5 0 
cm -3 and T kin = 60 — 120 K. The cross in Figure [ 6 ] 
marks the location of a representative good fit to all 
lines. Its values are n# 2 = 10 4,75 cm -3 and T k i n = 80 
K. T kin below 60 K challenges the observed ratios for 
H 13 CO+/H 13 CN and HN 13 C/H 13 CN , and T kin > 120 K 
begins to disagree with HCO + /H 13 CO + . Strictly speak¬ 
ing the HCO+/HCN ratio is not matched, however the 
disagreement is small (^30 %), and at the level that we 
expect other complications such as line width variation 
and transition depend filling factors to become impor¬ 
tant. The degree of internal consistency obtained in these 
simple uniform physical condition models is encouraging. 

The favored densities and kinetic temperatures of the 
dense gas isotopologues match well those previously 
found in Knudsen et al. (2007), with a range overlap¬ 


ping but favoring about 0.25 dex lower densities. These 
results confirm Knudsen et al. (2007)’s conclusions that 
HCN, HNC and (to a somewhat lesser degree) HCO + 
are moderately subthermally excited. The dense gas iso¬ 
topologues are even more strongly subthermal. For the 
above favored model, we obtain T ex = 34 K and 44 K 
for HCN and HCO + , respectively. Corresponding opac¬ 
ities r are 13 and 2.4 for HCN and HCO + . For the 
isotopologues, T ex [r] are 7.2 K [2.4], 16 K [0.27] and 5.3 
K [0.71] for H 13 CN, H 13 CO + and HN 13 C, respectively. 
LVG modeling demonstrates that while the main tran¬ 
sition line opacities crudely match those predicted from 
LTE, the isotopologue opacities are significantly larger 
than 12 r/ A. This stems from the fact that the effective 
critical density of a transition is dependent on its opac¬ 
ity (radiative trapping), so the isotopologues are more 
strongly subthermal. As a result, the level populations 
of the isotopologues settle into the lowest J transitions 
more efficiently, raising their opacity. This explains the 
well known effect that isotopic line ratios are gas den¬ 
sity probes (e.g., Meier et al. 2001). Furthermore this 
is strong evidence that the main species, particularly of 
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Fig. 7.— Same as Fig. [ 5 ] but for PDR tracer line ratios (Sec. |4.3| >. 

HCN, must have opacities well in excess of unity. 

The slightly larger HCN/H 13 CN and HCO+/H 13 CO+ 
ratios together with somewhat lower HN 13 C/H 13 CN ra¬ 
tios toward the outer nuclear disk are consistent with 
warm gas that is ^0.5 dex lower in density. However, 
the decrease in the H 13 CO + /H 13 CN ratio toward the 
outer disk is somewhat unexpected as gas density is de¬ 
creasing in this region. Raising T kin significantly will 
lower the H 13 CO + /H 13 CN ratio slightly even at lower 
density but only by ^10 % for T kin up 300 K. This 
effect does not appear to be strong enough to explain 
the observed ratio, and so this ratio is evidence for a 
small, but genuine abundance enhancement of HCN and 
its isotopologue relative to HCO + and HNC across the 
nuclear disk — or conversely a decrease of both the iso- 
topologues of HCO + and HNC relative to H 13 CN. In¬ 
voking mechanical heating may be a viable candidate for 
elevating HCN at the expense of HCO + and HNC (Loe- 
nen et al. 2008, Mejierink et al. 2011, Rosenberg et al. 
2014), because hotter gas converts HNC to HCN (e.g., 
Schilke et al 1992). If mechanical heating is more pro¬ 
nounced at the base of the outflow in the outer disk then 
this may explain the further drop in H 13 CO + /H 13 CN 
here. However this possibility is speculative so we con¬ 
sider the explanation for the lower H 13 CO + /H 13 CN ratio 
unsettled. 

4.3. PDR Tracers 

Given the strong radiation field in NGC253’s center, 
the presence of PDR tracers are expected and observed. 
PDRs warm gas and maintain significant amounts of car¬ 
bon in ionized form (e.g., Tielens & Hollenbach 1985), so 
tracers include species that form rapidly from C + and 
H or H2. These include simple hydrocarbons, such as 
CH, C2H, C-C3H2 and molecules that form directly from 
these species, like CN (e.g. Fuente et al. 1993, Sternberg 


& Dalgarno 1995, Rollig et al. 2007). In molecular gas 
strongly irradiated by UV photons we expect such species 
to exhibit elevated abundances. The two main PDR trac¬ 
ers in our dataset are CN and C2H. C2H has a critical 
density that is significantly lower than CN (see Table 4), 
and therefore traces somewhat more diffuse PDRs. 

Both CN and C2H exhibit fine/hyperfine structure 
splitting and so their opacity can be determined directly. 
For C 2 H the theoretical LTE C 2 H(3/2-l/2)/C 2 H(l/2- 
1/2) ratio should be 2.3 in the optically thin limit (for 
the subset of blended lines included), dropping to 1 
when optically thick. We observe a global ratio of ^2.4 
(Fig. [7]). Likewise for CN, the CN(3/2-l/2)/CN(l/2-l/2) 
fine structure ratio is 2 (Fig. [7]), which is what we observe 
for the inner nuclear disk of NCC 253. Both CN and C2H 
are likely optically thin and therefore their intensities are 
proportional to their column densities. 

The maps of these two species are strongly dominated 
by the innermost disk locations, 4-7, where star forma¬ 
tion is most intense and where the ionized/molecular 
outflow originates. Fig. [7] displays both the CN(l-0; 
3/2-l/2)/CO(l-0) (hereafter CN/CO) and the C 2 H(1- 
0;3/2—1/2)/HCN(l-O) (hereafter C 2 H/HCN) line ra¬ 
tio^] Both ratios decrease significantly from the in¬ 
ner (CN/CO — 0.11; C2H/HCN ~ 0.45) to the outer 
(CN/CO - 0.035; C 2 H/HCN ~ 0.18) nuclear disk. This 
behavior is explained by a combination of an increased 
PDR fraction and differential optical depth. To sep¬ 
arate these effects, ratios were calculated between op¬ 
tically thin proxies for CO and HCN, CN/C 17 0 and 
C 2 H/H 13 CN (Fig. 0 . In both cases the trend for en¬ 
hanced PDR tracers in the inner disk is seen, though 

11 These two ratios are chosen instead of the more common ratio, 
CN/HCN (e.g., Boger & Sternberg 2005), because only they can 
be obtained from matched array configurations in our dataset. 
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less pronounced. The average CN/C 17 0 ratio is <20 in 
the outer disk and a factor of ^2 higher towards the in¬ 
ner disk region). Likewise, for C2H/H 13 CN, the outer 
disk average is ^2 and the inner disk average is 3.5- 
4.5. The fact that C2H/H 13 CN and CN/C 17 0 both show 
the same magnitude of an effect suggests that density 
changes across the nucleus are not controlling the line 
ratio. Hence it appears that PDRs constitute a larger 
fraction of the dense gas towards the inner nuclear disk. 
Simple hydrocarbons exhibit elevated abundances in this 
environment. We note that the more complex hydrocar¬ 
bon symmetric top, CH3C2H, has a rather similar mor¬ 
phology to CN and C2H, suggesting that its abundance 
can also be elevated in PDRs. 

We also note that the PDR tracers are observed in the 
base of the molecular outflow. This is consistent with 
the recent detection of C2H in the molecular outflow in 
the nucleus of Maffei 2 (Meier & Turner 2012). 

4.4. Shock Tracers 

Given the intense star forming environment in 
NGC253’s center and its associated outflow of molec¬ 
ular gas (B13), we expect the presence of shocked gas 
tracers. Our dataset includes two common extragalac- 
tic shock tracers, SiO (e.g., Garcia-Burillo et al. 2000) 
and HNCO (Meier & Turner 2005). HNCO is thought 
to come from dust ice mantle sublimation (where due 
to filling factor reasons ice mantle evaporation in hot 
cores is invisible on larger scales), while SiO is enhanced 
by ejection of significant Si from sputtered silicate grain 
cores in high velocity shocks (e.g., Meier & Turner 2005, 
Rodriguez-Ferndandez et al. 2010, Tideswell et al. 2010, 
Mart in-Pint ado et al. 1992). 

SiO emission is distributed across the nucleus in a se¬ 
ries of compact sources with some relative enhancement 
seen toward the eastern outer nuclear disk. Calculated 
SiO abundances (X(SiO) ~ 10 -9 ) confirm that SiO is 
elevated across the entire nuclear disk. This SiO value is 
comparable to that seen locally on GMC scales in ‘shock 
regions’ in other nearby spiral galaxies (e.g., Usero et 
al. 2006, Meier & Turner 2012) and slightly larger than 
found in previous NGC 253 SiO measurements (differ¬ 
ences are due to different adopted T ex and N(H2) values; 
Garcia-Burillo et al. 2000). These abundances are much 
larger than observed toward quiescent Galactic molecular 
clouds (< 10 —11 ; e.g., Zuirys et al. 1989, Martin-Pintado 
et al. 1992). 

Interestingly, HNCO has the most distinctive morphol¬ 
ogy of all the bright lines, being completely dominated 
by the outer locations of the disk. Such an enhance¬ 
ment at the outer edges of the nuclear disk is also seen 
in several other nearby barred nuclei, IC342, Maffei 2 
and NGC 6946 (Meier & Turner 2005; 2012). Estimated 
HNCO abundances reach ~ 10 -8 across the outer nuclear 
disk and decrease to ^ 10 -9 towards the inner nuclear 
disk. HNCO abundance in quiescent Galactic clouds ap¬ 
proach this inner disk value (e.g. Marcelino et al. 2009, 
Tideswell et al. 2010). 

Assuming SiO traces shocks, then the SiO/HCN line 
ratio can be considered to first order an indicator of 
the fraction of the dense gas that is experiencing strong 
shocks. HNCO/SiO ratio is often considered a tracer of 
shock strength, because the shocks need to be stronger 
to elevate the SiO abundance (e.g., Usero et al. 2006). 


In Fig. [8] we present these two ratios across the nuclear 
disk of NGC 253. The SiO/HCN intensity ratio is fairly 
flat across the disk with slight elevations seen towards re¬ 
gions 2 and 3. This indicates that shocks are uniformly 
present over much of the nuclear disk in NGC 253. The 
HNCO/SiO ratio in the outer nuclear disk matches the 
ratio observed in bar shock regions of nearby spirals like 
IC 342 (Meier & Turner 2005, Usero et al. 2006) and Maf¬ 
fei 2 (Meier & Turner 2012) well. Towards the inner disk 
the ratio drops dramatically to a value of ^0.17. Only 
M 82 has a ratio anywhere near this low (Garcia-Burillo 
et al. 2000, Martin et al. 2009). 

At face value this extremely low HNCO/SiO ratio in 
the inner disk could be taken to indicate that shock 
strengths decrease away from the core. However, the SiO 
morphology and abundance suggests that, even in the 
outer disk, shock strengths are strong enough to keep 
SiO high relative to HCN. Therefore other effects that 
further alter the HNCO morphology besides simple shock 
strength must be considered. Two important second or¬ 
der effects that can alter the HNCO/SiO ratio (see Meier 
& Turner 2012 for a thorough discussion): 1) HNCO has 
a significantly higher photodissociation rate than SiO, so 
its shock signatures are erased more easily in the pres¬ 
ence of PDRs/UV radiation fields. 2) SiO is a linear top 
molecule so its partition function is oc T ex , while HNCO 
is an asymmetric top molecule, with a partition func¬ 
tion oc T^ 5 . So in hot gas, HNCO’s lower energy states 
depopulate more rapidly than SiO. For HNCO this ef¬ 
fect is mitigated at T ex ~ 5 — 20 K by the fact that 
E w (HNCO) is ~ 11 K compared to SiO’s value of 6.3 K, 
but if T ex ~ 75 K in the inner disk then this could be a 
pronounced effect. 

In the inner n ucle ar disk, the gas is hot and PDR dom¬ 
inated (see Sec. 4.3), so we expect that both act to drive 
down the HNCO/oiO ratio in the center. Therefore it 
appears we are observing the starburst (and outflow) ac¬ 
tively in the process of erasing the signatures grain man¬ 
tles species in the inner disk (Martin et al. 2009a,b). In 
this context it is not surprising that its HNCO/SiO ratio 
is similar to M 82’s, a nucleus where PDRs have com¬ 
pleted the task of erasing the weak shock tracers (Takano 
et al. 1995, Garcia-Burillo et al. 2001) 


4.5. Sulfur Species 

We detect a significant number of sulfureted species, 
including CS, SO, NS, H 2 CS, CH 3 SH and C 2 S. When 
accounting for the different array configurations used to 
observe them, CS, SO, H 2 CS and NS all have quite sim¬ 
ilar morphologies. They tend to be dominated by the 
dense GMCs, with a preference toward the clouds on the 
northeast side of the continuum peak (7 and 8), whereas 
the PDR/hydrocarbons tend to favor the southwestern 
side of the inner disk (position 5). Interestingly, the mor¬ 
phology of the sulfuretted species, particularly CS and 
SO are v ery s imilar to SiO. Assuming SiO traces shocked 
gas (Sec. |4.4| ), this hints at a possible connection between 
the S species and shocked gas. This is somewhat unex¬ 
pected, at least for CS, given that in other nearby nuclei 
it has been argued that CS is primarily a PDR tracer 
(e.g., Meier & Turner 2005; Martin et al. 2009). Ev¬ 
idently a mix of shocks and strong radiation fields are 
capable of maintaining high abundances of SiO and CS, 
though not HNCO. 
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Fig. 8.— Same as Fig. [ 5 ] but for shock tracer ratios (Sec. |4.4| ) 
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Fig. 9. — Same as Fig. [5] but for sulfur species line ratios (Sec. |4.5] >. 

The CS/SO abundance ratio is sensitive to the C + /0 
abundance ratio and it has been argued to be either a 
tracer of the C/O abundance ratio (~ 0.4 in the Galaxy), 
or early vs. late time chemistry (e.g., Heithausen et al. 
1998, Nilsson et al. 2000). So when C or C + is nearly 
as abundant as atomic O, CS is strongly favored over 
SO. Our observed SO/CS intensity ratio is ^0.1-0.18 
over much of the nuclear disk (Fig. [9]). This corresponds 
to a CS/SO abundance ratio of >L5-3 (higher if CS is 
optically thick), in good agreement with Galactic star 
forming regions. Nilsson et al. (2000) find values ranging 
between 0.4-5, with most at ^2. 

Likewise, the NS/SO (and NS/CS) abundance ratios 
are particularly sensitive to atomic C and O abundances. 

NS is formed from the neutral-neutral reactions, S + NH 
—» NS + H and N + SH NS + H. Atomic O destroys 
both NS and its parent SH, converting them to NO and 
SO, respectively. Likewise, reactions with C/C + drive 
NS into CS and CN. In star forming regions the radia¬ 
tion field can keep a significant amount of C and O in 
non-molecular form. As such, standard gas-phase models 
predict NS cannot not attain abundances above X(NS) 

~ 10 -11 (e.g., McGonagle & Irvine 1997, McElroy et al. 
2013). In the inner nuclear disk we observe NS abun¬ 
dances in the range of 1 - 6xl0 -9 , at least two orders of 
magnitude larger than the maximum gas-phase models 
can accommodate. NS/SO abundances ratios are also 
high, reaching 0.2 - 0.5 (the exact values are uncertain 


due to non-matched array configurations). These ele¬ 
vated abundances are consistent with that determined 
from single-dish measurements toward NGC 253 (Martin 
et al. 2003, 2005), except the present observations further 
demonstrate that the high abundances are attained to¬ 
ward the same inner disk GMCs where PDR conditions 
are pronounced. 

To explain the similar (high) abundances of CS, SO 
and NS toward the same region requires non-standard 
chemistry. The atomic abundances of species like O must 
be kept very low, despite the PDR conditions. Two pos¬ 
sible scenarios appear feasible. The first, introduced by 
Viti et al. (2001) (see also Charnley 1997, Harada et al. 
2010 and Izumi et al. 2013), invokes a high tempera¬ 
ture gas phase chemistry due to the passage of a shock. 
This high temperature gas drives O into the more inert 
H 2 0 form so it does not destroy NS as rapidly. A sec¬ 
ond possibility is that in the pre-starburst phase of NGC 
253’s nucleus, much of the atomic species condensed onto 
the dust grains where they are converted into their sat¬ 
urated form (e.g., Brown et al. 1988). At a later time 
shocks liberated the saturated species (e.g., CH4, NH3, 
H 2 0 and possibly H 2 S). The gas-phase chemistry then 
evolved from molecular initial conditions, keeping atomic 
C and O low. Both of the scenarios invoke the presence 
of shocks, and so are consistent with the tight morpho¬ 
logical connection between the sulfur species and SiO. 

We observe both ortho-H 2 CS at 101 GHz and a mix 
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of para transitions at 103 GHz. The ortho-to-para ra¬ 
tio of H 2 CS, the sulfur analog to formaldehyde H 2 CO, is 
useful for constraining the formation conditions of H 2 CS 
(Minh et al. 1991). Under optically thin, high-T (>15 K 
in the case of H 2 CS), LTE case, the ratio of the ortho 
species of H 2 CS (those with K a [K_i] even, see Tab. fl]) 
to the para (those with K a odd) should be 3.0. If the 
temperature at the time of formation is <15 K then the 
two forms of H 2 CS with not be thermalized to their sta¬ 
tistical weights and hence have an abundance ratio less 
than 3 — approaching unity at T=0K. The intensity 
ratio of the two types of transitions (Fig. |9| is approxi¬ 
mately constant across the disc at a value of ^1.4. For 
T e £ = 75 K the measured ratio implies an ortho-to-para 
ratio of 2.9. Therefore every position is consistent with 
standard LTE, high temperature limit. Thus, at least 
in the case of H 2 CS, this suggests that H 2 CS formed in 
the high temperature gas phase after evaporation, or if 
it was formed on the grain surfaces in an earlier epoch, 
those surfaces were still significantly warmer than 15 K. 

4.6. Millimeter Recombination Lines 

We have detected three hydrogen recombination lines 
in our observations (H40cq H50/? and H52/?). The mor¬ 
phology of these recombination lines is very different 
from all the molecular gas tracers (e.g., Figs. IT] and p|: 
The emission is centrally concentrated, coincident with 
the location of the underlying continuum emission (L14). 
H (3 radio recombination lines are reported outside the 
Local Group for the first time here. The error bars are 
significant for the (3 lines, but the global H40<a/H50/3 ra¬ 
tio is around 3.5, in good agreement with the LTE value 
of 3.5 (=1/0.285) for n e = 10 4 cm' 3 and T e = 10 4 K. 
The 7 lines would be fainter than the /? line by a factor of 
2.3. Indeed, we do not clearly detect any 7 lines. We also 
do not detect any He<a lines. Other radio/mm recombi¬ 
nation lines mapped in NGC 253 have been discussed ex¬ 
tensively elsewhere (e.g., Anantharamaiah & Goss 1996, 
Mohan et al. 2002, 2005, Rodriguez-Rico et al. 2006, 
Kepley et al. 2011). The morphology of both the a and 
/? lines mapped here are in excellent agreement with the 
ones reported previously. 

4.7. Tentatively Identified Species 

We here briefly comment on other selected species de¬ 
tected in our observations. 

A large number of tentatively identified lines can 
be attributed to carbon-rich complex species (e.g., 

ch 3 ch 2 cn, ch 3 ch 2 oh, ch 2 cn, ch 2 co, ch 3 cho 

and NH 2 CHO (Tab. [l]). These species represent the 
next step in complexity beyond the molecules detected 
in external galaxies up until now, and suggest that these 
molecules may be fairly extended in NGC 253. 

One unidentified feature detected is of particular inter¬ 
est. This is the feature U-114.22 at ^114.218 GHz. This 
spectral feature is within <4 MHz (<11 km s -1 ) of both 
C 4 H(N=12-11; J=23/2-21/2; F=12-ll) and the v = 1 
transition of CO(1-0). The upper energy state of the 
CO(1-0; v = 1) line is 3089 K above ground, so the exci¬ 
tation temperature of the gas must be >450 K in order to 
explain the observed line ratio. The morphology of this 
transition is very similar to the RRLs (accounting for 
resolution differences), as would be expected if it is trac¬ 
ing molecular gas not associated with the cold molecular 


phase. The C4H transition has a much lower upper en¬ 
ergy level of 35.6 K, typical of the other transitions seen 
here. Moreover, it is plausible to expect that the C4H 
might trace PDR emission (like C 2 H), which are also con¬ 
centrated towards the inner disk, though why it would be 
significantly more compact than C 2 H is unexplained. We 
consider C4H as the most likely identification, however 
the observed morphology indicates that the possibility 
this feature could be vibrationally excited CO deserves 
further investigation. 

5. SUMMARY AND CONCLUSIONS 


We present the detection of a total of 50 molecular lines 
in the nearby starburst galaxy NGC 253, based on early 
science (cycle 0) observations taken with the Atacama 
Large (Sub-)Millimeter Array (ALMA). For 27 lines we 
have an unambiguous identification, and 13 lines have 
tentative IDs (no plausible line identification was possible 
in 10 cases). 

We here describe a schematic picture of the central 
starburst region of NGC 253 mapped by ALMA (Fig. 10): 
As described in Sec. |3.1[ the nucleus of NGC 253 is cnar- 
acterized by the inner portion of NGC253’s large-scale 
bar. The outer part of the nuclear disk (the ‘outer nu¬ 
clear disk’), likely represents the location where the gas 
is flowing radially inward along the large-scale bar. In¬ 
side this is a compact region exhibiting a large quantity of 
high density gas and intense star fo rmat ion, as evide nced 
by the presence of dens e ga s (Sec. 4.2), PDR (Sec. 4.3) 
and shock tracers (Sec. |4 4| ). This component (the mi¬ 
ner nuclear disk’), dominates the morphology of most of 
the spectral lines detected by the observations. We de¬ 
tect high density molecular gas tracers (HCN, HCO + and 
CN) at the base of the molecular outflow first detected 
in the 12 CO emission (B13). 

In detail, we find moderate 12 CO opacities of ~2—5, 
despite the large column density towards NGC253’s nu¬ 
cleus. This may be due to increased line-widths due 
to turbulence, non-circular motions and/or elevated gas 
temperatures that reduce the CO opacity per unit veloc¬ 
ity. Comparing HCN(l-O) and HCO + (1-0) with their 
13 C-substituted isotopologues yields that the HCN(1- 
0) and HCO + (1-0) are also optically thick, with similar 
(moderate) opacities to CO. These high HCN and HCO + 
opacities imply that the main isotopologue HCN/HCO+ 
(and HCN/HNC) line ratios have less diagnostic power in 
this starburst environment. Using the isotopically sub¬ 
stituted versions of these high density tracers we measure 
much more robust values for the popular ”HCO + /HCN”, 
’’HNC/HCN” and ”HNC/HCO + ” ratios. All five of the 


isotopic dense gas tracer ratios can be well fit by LVG 
models with uh 2 = 10 4,75 cm -3 , T kin = 80 K and 
X(HCN) ~ 5 X(HCO + ) ~ 5X(HNC). 

The weak shock tracer/ice mantle species HNCO has 
the most distinctive morphology of all the bright lines be¬ 
ing completely dominated by the outer locations of the 
disk. This is at odds with what is seen in the strong 
shock tracer, SiO. Attributing the dramatic variation in 
the HNCO and SiO maps across the nucleus to chang¬ 
ing shock strength is unsatisfactory because SiO is en¬ 
hanced across the nucleus at a level that should imply 
shocks are strong throughout the inner disk. A possible 
explanation for the faintness of HNCO is that because 
of its higher photodissociation rate relative to SiO, the 
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Fig. 10.— A schematic picture of the central starburst region of NGC 253 (see discussion in Sec. [5j . The morphology is based on Figure 
1 of (B13) with the chemistry added. 


ice mantle evaporation shock signatures are being pref¬ 
erentially erased in the presence of the intense central 
radiation fields. This is consistent with the presence of 
an enhanced PDR fraction towards the region of HNCO’s 
dramatic decrease as determined form the CN/C 17 0 and 
C 2 H/H 13 CN ratios. The sulfureted species appear also 
to be (indirectly) connected to the presence of shocks 
across the inner disk. 

We have detected three hydrogen recombination lines 
(H40<a, H50/3 and H52/?) that show a centrally concen¬ 
trated morphology (similar to the underlying continuum 
emission) that is distinctly different from all molecular 
gas tracers. 

Finally, we report some tentative identifications, in¬ 
cluding the aldehydes CH 3 CHO and NH 2 CHO, cyanides 
CH 2 CN and CH 3 CH 2 CN and a number of organics of 
similar complexity, CH 2 CO, HCOOH and CH 3 CH 2 OH. 
The first vibrational state of the 12 CO(1-0) line is also 
possibly detected. If follow-up, multi-line studies con¬ 
firm these IDs then they would represent the first extra- 


galactic detections of the above species. 
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TABLE 3 
Line Intensities 


Transition 

Beam 

(") 

1(1) 

1(2) 

1(3) 

1(4) 

1(5) 

(K km 

1(6) 

s- 1 ) 

1(7) 

1(8) 

1(9) 

(110) 

SO(j JV =2 2 -l 1 ) 

2 . 0 / 2.4 

<5.7 

< 5.7 

< 7.2 

6.0 ± 1.8 

12.0 ± 3.6 

< 7.2 

19.5 ± 1.9 

< 5.7 

< 5.7 

< 5.7 

H 13 CN(i-o) 

2 . 0 / 2.4 

<5.7 

8.5 ± 2.6 

25.3 ± 2.5 

39.6 ± 4.0 

58.4 ± 5.8 

< 7.1 

64.3 ± 6.4 

45.3 ± 4.5 

11.7 ± 1.2 

13.4 ± 1.3 

HCO(j k<j k c =1 i ,o _ o 0 ,o) 

1 . 9 / 2.4 

<5.6 

< 5.6 

< 7.1 

8.2 ± 2.5 

< 7.1 

< 7.1 

< 7.1 

< 5.6 

< 5.6 

< 5.6 

H 13 CO+(i-o) 

1 . 9 / 2.4 

<5.6 

< 5.6 

8.7 ± 2.6 

37.8 ± 3.8 

41.7 ± 4.2 

11.0 ± 3.3 

49.1 ± 4.9 

26.8 ± 2.7 

< 5.6 

7.3 ± 2.2 

SiO ( 2 —1; v=o) 

1 . 9 / 2.4 

<5.6 

15.8 ± 1.6 

36.6 ± 3.7 

41.6 ± 4.2 

48.1 ± 4.8 

< 7.0 

62.7 ± 6.3 

46.4 ± 4.6 

13.2 ± 1.3 

14.3 ± 1.4 

HN 13 C(i-o) 

1 . 9 / 2.4 

<5.6 

5.9 ± 1.8 

< 7.0 

17.3 ± 1.7 

22.1 ± 2.2 

7.3 ± 2.2 

26.6 ± 2.7 

13.1 ± 1.3 

< 5.6 

< 5.6 

C2H(at = 1-0; J=3/2-l/2) 

1 . 9 / 2.4 

<5.6 

27.3 ± 2.7 

46.2 ± 4.6 

194.8 ± 19.5 

224 ± 22 

157± 16 

246 ± 25 

133 ± 13 

36.5 ± 3.6 

30.9 ± 3.1 

C2H(l-0;l/2-l/2) 

1 . 9 / 2.4 

6.7T2.0 

15.0 ± 1.5 

19.4 ± 1.9 

67.2 ± 6.7 

95.3 ± 9.5 

69.7 ± 7.0 

98.9 ± 9.9 

49.7 ± 5.0 

14.4 ± 1.4 

12.8 ± 1.3 

HNCO(4 0 , 4 - 3 0 , 3 ) 

1 . 8 / 2.4 

<5.5 

43.1 ± 4.3 

67.8 ± 6.8 

25.5 ± 2.6 

7.1 ± 2.1 

< 5.5 

19.5 ± 1.9 

46.4 ± 4.6 

33.5 ± 3.3 

50.0 ± 5.0 

H(52)/3 

1 . 8 / 2.4 

<5.4 

< 4.1 

< 5.4 

< 4.1 

9.0 ± 2.7 

14.9 ± 1.5 

< 5.4 

< 5.4 

< 4.1 

< 4.1 

HCN(i-o) 

1 . 8 / 2.4 

24.0±2.4 

119 ± 12 

265 ± 27 

435 ± 44 

582 ± 58 

360 ± 36 

526 ± 53 

512 ± 51 

211 ± 21 

179 ± 18 

HCO+0-o) 

1 . 8 / 2.4 

26.7T2.7 

75.5 ± 7.5 

137 ± 14 

432± 43 

512 ± 51 

276 ± 27.6 

450 ± 45 

422 ± 42. 

176 ± 18 

153 ± 15 

CS(2-l) 

1 . 6 / 2.4 

<4.4 

57.2 ± 5.7 

99.7 ± 10.0 

208 ± 21 

233 ± 23 

136 ± 14 

287 ± 29 

225 ± 23 

63.4 ± 6.3 

53.0 ± 5.3 

H(40)a 

1 . 6 / 2.4 

<3.2 

< 3.2 

< 4.3 

< 3.2 

25.1 ± 2.5 

47.6 ± 4.8 

13.0 ± 1.3 

< 3.2 

< 3.2 

< 3.2 

H(50)/3 

1 . 6 / 2.4 

<4.3 

< 3.2 

< 4.3 

< 3.2 

7.8 ± 2.3 

12.1 ± 1.2 

< 4.3 

< 4.3 

< 3.2 

< 3.2 

SO(3 2 - 2,) 

1 . 6 / 2.4 

<4.3 

10.3 ± 1.0 

24.5 ± 2.5 

30.1 ± 3.0 

30.3 ± 3.0 

12.1 ± 1.2 

46.9 ± 4.7 

35.5 ± 3.6 

10.8 ± 1.1 

17.2 ± 1.7 

HC 3 N(ii-io) 

4 . 1 / 4.6 

<2.0 

16.2 ± 1.6 

29.4 ± 2.9 

60.7 ± 6.1 

59.9 ± 6.0 

52.8 ± 5.3 

72.0 ± 7.2 

48.0 ± 4.8 

23.2 ± 2.3 

19.7 ± 2.0 

CH 3 SH(4 0 - 3 0 )A+E 

4.2/4.6 

<2.0 

< 2.0 

< 2.3 

2.3 ± 0.7 

< 2.3 

< 2.3 

3.4 ± 1.0 

2.8 ± 0.8 

< 2.0 

< 2.0 

H 2 CS(3 1j 3 -2 1i2 ) 

4.2/4.6 

<2.0 

< 2.0 

3.1 ± 0.9 

4.1 ± 0.4 

3.1 ± 0.9 

3.4 ± 1.0 

4.8 ± 0.5 

4.8 ± 0.5 

2.6 ± 0.8 

2.8 ± 0.8 

CH 3 C 2 H(6 fe - 5 fc ) 

4.2/4.6 

<1.9 

5.8 ± 0.6 

5.4 ± 0.5 

24.2 ± 2.4 

24.8 ± 2.5 

20.9 ± 2.1 

21.5 ± 2.1 

11.9 ± 1.2 

6.2 ± 0.6 

4.4 ± 0.4 

H 2 CS 

4.1/4.6 

<1.9 

< 1.9 

< 2.2 

2.5 ± 0.8 

2.2 ± 0.7 

< 2.2 

3.8 ± 1.1 

2.8 ± 0.8 

2.2 ± 0.7 

2.6 ± 0.8 

C 17 0(i-o) 

4.0/4.6 

<2.7 

< 2.7 

3.9 ± 1.2 

9.2 ± 0.9 

10.2 ± 1.0 

9.5 ± 0.9 

9.1 ± 0.9 

7.0 ± 0.7 

5.0 ± 1.5 

4.2 ± 1.2 

CN( 1 - 0 ; 1 / 2 - 1 / 2 ) 

4.0/4.6 

4.9T1.5 

18.0 ± 1.8 

23.1 ± 2.3 

169 ± 17 

215 ± 22 

201 ± 20 

196 ± 20 

107 ± 11 

42.8 ± 4.3 

30.4 ± 3.0 

CN(l-0;3/2-l/2) 

4.0/4.6 

27.6T2.8 

58.0 ± 5.8 

52.9 ± 5.3 

322 ± 32 

403 ± 40.3 

380 ± 38 

325 ± 32 

206 ± 21 

109 ± 11 

81.9 ± 8.2 

CO(i-o) 

3.7/4.6 

718±72 

1920 ± 190 

1570 ±160 

3200 ± 320 

3550 ± 360 

3410 ± 340 

3130 ± 310 

2710 ± 270 

2160 ± 220 

1920 ± 190 

NS(5/2-3/2;7/2-5/2) 

3.7/4.6 

<3.2 

< 3.2 

< 3.9 

6.9 ± 0.7 

7.8 ± 0.8 

8.1 ± 0.8 

9.2 ± 0.9 

3.4 ± 1.0 

< 3.2 

< 3.5 


Note. — Integrated intensity measurements at positions 1—10 and intensity of the entire galaxy. Column 1: The transition, column 2: The original resolution of the ALMA data/common 
resolution used in this study, column 3-12: The integrated intensity for regions #1—#10 (see Fig. |4]for location of regions) with uncertainty. Uncertainties are set to 10%, representing absolute 
flux calibration uncertainties only, unless S/N is low then they are conservatively set to 30 % (see section [2j. Small errors associated with missing flux are not included. We caution that values 
for region 6 may be affected by absorption towards the central continuum source. All upper limits are 5cr. 
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TABLE 4 

LTE ’Reference’ Molecular Abundances for NGC 253 a 


Mol. 

Tier 0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

h 13 cn 

6.48 

<0.2-<l(-9) 

l- 8 (- 10 ) 

0.4-2(-9) 

0.4-2(-9) 

0.5-3(-9) 

<0.6-<3(-10) 

0.6-3(-9) 

cT 

1 

co" 

1 

10 

6 

2-9(-10) 

0.2-1 (-9) 

h 13 co+ 

5.33 

< 0 . 1 -< 8 (- 10 ) 

0.5-3(-9) 

0.8-4(-10) 

0 . 2-1 (-9) 

0.2-l(-9) 

0.6-3(-10) 

0.3-2 (-9) 

2-9(-10) 

<0.7-<4(-10) 

0.7-4(-10) 

SiO 

5.48 

<0.3-<l(-9) 

0.3-1 (-9) 

0.7-3.0(-9) 

0.5-2(-9) 

0.5-2(-9) 

<0.8-<3(-10) 

0.7-3(-9) 

0.6-3 (-9) 

0.2-1 (-9) 

0/3-1 (-9) 

hn 13 c 

5.62 

<0.3-<2(-9) 

l- 6 (- 10 ) 

<l-<7(-10) 

0.2-1 (-9) 

0.8-4(-9) 

0.8-4(-10) 

0.3-1 (-9) 

0.2-9(-10) 

<0.9-<5(-10) 

<<1-5 (-10) 

c 2 h 

5.27 

<0.4—<2(-8) 

1 —5(-8) 

2 - 8 (- 8 ) 

0.4-2(-7) 

0.4-2(-7) 

0.3-2(-7) 

0.5-3(-7) 

0.3-2(-7) 

l- 6 (- 8 ) 

1-5 (- 8 ) 

HNCO 

5.20 

< 0 . 1 -<l(- 8 ) 

0.4-3 (- 8 ) 

0.6-5(-8) 

0 . 1-1 (- 8 ) 

0.3-3(-9) 

<0.3-<2(-9) 

l-9(-9) 

0.3-2(-8) 

0.3-2(-8) 

0.5-4(-8) 

HCN 

6.42 

l-5(-9) 

0 . 2 -l(- 8 ) 

0.4-2(-8) 

0.4-2(-8) 

0.5-3(-8) 

0.3-2 (- 8 ) 

0.5-3(-8) 

0.6-3(-8) 

0.3-2(-8) 

0.3-2(-8) 

HCO+ 

5.37 

0.7-4(-9) 

0.7-4 (-9) 

l-7(-9) 

0 . 2-1 (- 8 ) 

0.3-1 (- 8 ) 

1-7 (-9) 

0.3-1 (- 8 ) 

0.4-2(-8) 

0.2-8 (-9) 

l-7(-9) 

CS 

5.78 

<0.5-<2(-9) 

2-9(-9) 

0.4-2(-8) 

0.5-2(-8) 

0.5-2(-8) 

0.3-l(-8) 

0.7-3(-8) 

0.6-3(-8) 

0.3-l(-8) 

2-9 (-9) 

so 

5.47 

<2-<7(-9) 

2-6(-9) 

0.4—2(-8) 

0.3-1 (- 8 ) 

0.3-1 (- 8 ) 

1-4 (-9) 

0.5-2 (- 8 ) 

0.4-2(-8) 

2-6(-9) 

0.3-1 (- 8 ) 

hc 3 n 

5.71 

<l-<0.7(-10) 

3-2(-10) 

6-4(-10) 

7-5(-10) 

7-4(-10) 

6-4 (-10) 

9-5(-9) 

7-4(-10) 

4-3(-10) 

4-2 (-10) 

CH 3 SH 


<l-<7(-9) 

<0.4-<3(-9) 

<0.5-<3(-9) 

0.3-2(-9) 

<0.2—<2(-9) 

<0.3-<2(-9) 

0.4-3 (-9) 

0.4-3(-9) 

<0.3-<2(-9) 

<0.4-<3(-9) 

h 2 cs 

5.25 

<2-<6(-9) 

<0.8-<2(-9) 

<0.9-<3(-9) 

0.6-2(-9) 

0.5-1 (-9) 

0.5-1 (-9) 

0.9-3(-9) 

0.8-2(-9) 

0.8-2 (-9) 

1-3 (-9) 

CN 

6.24 

0.4-2(-8) 

0.4-2 (- 8 ) 

0.4-2(-8) 

l- 6 (- 8 ) 

l- 6 (- 8 ) 

2-7(- 8 ) 

l-7(-8) 

l-5(-8) 

0.6-3 (- 8 ) 

0.5-2(-8) 

NS 


<0.3-<l(-8) 

<l-<4(-9) 

<l-<4(-9) 

l-5(-9) 

1—5(-9) 

2-5(-9) 

2-6(-9) 

0.8-3(-9) 

<l-<3(-9) 

<l-<4(-9) 

c 17 o 

3.28 

<0.6-<3(-7) 

<0.2-<l(-7) 

0.3-2(-7) 

0.4-2(-7) 

0.4-2(-7) 

0.4-2(-7) 

0.4-2(-7) 

0.4-2 (-7) 

0.4-2(-7) 

0.3-2(-7) 

NpD 


3.6(22) 

9.6(22) 

7.8(22) 

1.6(23) 

1.8(23) 

1.7(23) 

1.6(23) 

1.4(23) 

1.1(23) 

9.6(22) 


a Entries have the form: a(b) = a X 10 b . Molecular hydrogen column densities are determined from the CO(l-O) intensity assuming a conversion factor of Xco = 0.5 X 10 20 cm -2 (K km s -1 ) -1 
(Bolatto et al. 2013a). Calculation of the column densities of individual species are calculated assuming LTE excitation. The range quoted for each entry corresponds to 10 - 75 K. Overall 
systematic uncertainties, including the adopted excitation temperature (range shown), galactic position (positions 1-3 and 7-10 likely have excitation temperatures at the low end of the range), 
and the CO conversion factor (uncertain to a factor of three; Bolatto et al. 2013a), dominate the error budget, so separate error bars are not presented for each position. Upper limits are 5cr. 
b The critical density (not including radiative trapping effects), log(n cr = Aij/Cij [100 K]). Rates are adapted from the Leiden LAMDA database, van der Tak et al. 2007, with Cij’ s: 
HCN/H 13 N/HN 13 C: Dumouchel et al. (2010); HCO+/H 13 CO+: Flower (1999); SiO: Dayou & Balanca (2006); C 2 H: Spielfiedel et al. (2012); HNCO: Green (1986); CS: Lique et al. (2006a); SO: 
Lique et al. (2006b); HC 3 N Green & Chapman (1978); H 2 CS: Wiesenfeld & Faure (2013); CN: Lique et al. (2010); C 17 0: Yang et al. (2010). 
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